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Abstract 

A theoretical study is made of the magnetization versus applied field curves of ferro- 
magnetic/nonmagnetic multilayers constructed according to a Fibonacci quasiperi- 
odic sequence. The ferromagnetic films are assumed to have uniaxial anisotropy 
and are coupled by both bilinear and biquadratic effective exchange. The effects 
of quasiperiodicity in the magnetic phases are illustrated numerically for Fe/Cr 
systems. 
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Magnetic multilayers have attracted attention due to their physical properties 
and potential for technological applications [1]. The coupling between the 
magnetic layers often corresponds to a bilinear (BL), or Heisenberg, exchange 
type. However, some multilayers also exhibit a biquadratic (BQ) exchange 
coupling, which favours a 90° orientation between the adjacent ferromagnetic 
layers [2,3]. The BQ exchange, although often weaker than the BL exchange, 
can be large in some systems. This leads to new possibilities for the phase 
diagrams and magnetization curves. 

Furthermore, it is possible to build magnetic multilayers where the film thick- 
ness and stacking are well controlled, as in periodic superlattices. However, 
there has also been interest in magnetic systems formed with the layers fol- 
lowing a quasiperiodic sequence. Recent studies, applied to Fe/Cr Fibonacci 
multilayers, showed that a self-similar behaviour of the magnetic properties 
occurs only when the BQ and BL exchange are comparable [4,5]. The authors 
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included effects of a cubic anisotropy (see, e.g., [6]). However, some magnetic 
multilayers may instead exhibit uniaxial anisotropy [7], and this gives a new 
symmetry to the system that, together with the BQ exchange and quasiperi- 
odicity, may lead to new configurations. This has motivated the present work. 

A Fibonacci structure is formed by juxtaposing two basic building blocks, 
A and B, following a Fibonacci sequence. For the present work we choose 
ferromagnetic Fe and nonmagnetic Cr to correspond to A and B, respectively. 
By definition, the Nth Fibonacci generation, Sn, is obtained by appending 
generation N — 2 to generation N — 1, i.e., Sn = Sn-iSn-2 taking So = B and 
Si = A. Thus, for example, S3 and S$ correspond to the magnetic multilayers 
Fe/Cr/Fe and Fe/Cr/Fe/Fe/Cr/Fe/Cr/Fe respectively. 

We assume the magnetization of the ferromagnetic layers to be in the xy- 
plane and take the z-axis as the growth direction. Following [4,5] but with the 
anisotropy replaced by a uniaxial term as in [7] , the magnetic energy per unit 
area is 
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Index i(= 1,2, ... ,n) labels the Fe layers in the Nth generation structure. For 
example, n = 4 in the case of N = 5, with three Fe layers having thickness 
t and one having thickness It, where t denotes the thickness of a single A 
layer (the basic tile). In layer % (with thickness tj) the equilibrium magneti- 
zation is in the xy-plane at angle 0, to the x-axis. The first term in eq. (1) 
is the Zeeman energy due to an applied field H at an angle 9h- The second 
term is due to a uniaxial anisotropy H ua with easy axis specified by angle 6 ua . 
The third and fourth terms (involving parameters Hu, Hf, q ) describe the BL 
and BQ exchange between adjacent layers. The set {6} t of equilibrium an- 
gles is calculated numerically by minimizing the magnetic energy given by eq. 
(1). For that purpose we use two methods, namely, the simulated annealing 
and the gradient methods (see [5] for descriptions). The simulated annealing 
method is based on the fact that heating and then cooling a material slowly 
brings it into a more uniform state, which is the minimum energy state. In this 
process, the role played by a pseudo temperature T is to allow the configura- 
tions to reach higher energy states with probability p given by the Boltzmann 
law p = exp(—AE/kT), where AE is the energy difference. Energy barriers 
that would otherwise force the configurations into local minima can then be 
overcome. On the other hand, the gradient method is based on finding the 
directional derivative of the magnetic energy in the search for its global min- 



imum in the n-dimensional space composed of the variables {0} i - It is the 
gradient of the magnetic energy with relation to the angles that furnishing 
the direction, and eventually the location, of the required global minimum. 
Both cited methods are used for each value of the applied magnetic field and 
for each set of magnetic parameters. We choose the configuration with the 
lowest energy furnished by both methods as giving the equilibrium configu- 
ration {0} i . From the calculated {0} i , we obtain normalized values of the net 
magnetization in the H direction, appropriate to the Fe/Cr system. 

Magnetization curves found for N = 3,5 and 7 are shown in Fig. 1. Here 
the in-plane magnetic field is applied along the easy axis. The magnetic pa- 
rameters are such that Hb q << \Hm\ and Hb q << H ua . For N = 3 (see Fig. 
la) the magnetizations in adjacent layers are antiparallel at low fields. As H 
increases, a first-order phase transition occurs (at ~ 0.8 kOe) to a spin-flop 
phase. Saturation is reached at H ~ 1.9 kOe. For N = 5 (see Fig. lb), due to 
the different thickness of the Fe layers, the net magnetization at low fields has 
about 20% of its saturation value. There are first-order phase transitions at 
H ~ 1.0 kOe and H ~ 1.6 kOe. Saturation is reached at ~ 2.9 kOe. When 
N = 7 (Fig. lc) there is again a non-zero magnetization at low fields. There 
are three first-order phase transitions at H ~ 0.5 kOe, 1.0 kOe and 1.6 kOe. 

Another set of magnetization curves is shown in Fig. 2, taking a larger BQ 
exchange and other parameters as before. For N = 3, due to the strong Hb q the 
adjacent layer magnetizations are only approximately antiparallel at low field 
(corresponding to an asymmetric phase [7]). A first-order phase transition 
occurs at Hq ~ 0.4 kOe to a spin-flop phase and saturation is reached at 
Ho ~ 4.3 kOe. For N = 5 all transitions appear continuous (second-order 
phase transitions) and saturation is reached at Hq ~ 6.4 kOe. For N = 7 there 
appears to be a first-order phase transition at H ~ 0.15 kOe and saturation 
is reached at H ~ 6.6 kOe. Again the fifth and seventh generations exhibit a 
nonzero magnetization in the low-field region. 

We have also studied cases where the applied field does not lie along the easy 
axis (6h 7^ ua ), and an example for N = 3 is shown in Fig. 3. When the field 
is applied at 45° and 90° from the easy axis, the magnetizations continuously 
rotate towards the field direction. There is no first-order phase transition in 
these configurations. The saturation fields for 45° and 90° are H$ ~ 2.3 kOe 
and 2.9 kOe, respectively. 

In conclusion, we have modelled the magnetization versus applied field curves 
of Fe/Cr Fibonacci multilayers including both BL and BQ exchange. The role 
of uniaxial anisotropy is studied, including situations when the applied field 
direction is different from the easy axis. One can see a discrete self-similar 
pattern in the system when Hb q « \Hu\ (Fig. 1) and most transitions are 
of the first-order type. This property is highlighted by the boxes drawn as 



insets in Fig. 1. Specifically, the magnetization curves of higher generation 
N reproduce some aspects of the magnetization curves of lower generation 
N — 2. By contrast, when H^ ~ \Hm\ and most of the transitions are second- 
order type, there is no apparent self-similarity This behaviour, regarding the 
effect of BQ exchange, differs from that in quasiperiodic multilayers having 
cubic anisotropy [4,5]. We may surmise that the lower symmetry of uniaxial 
anisotropy changes the nature of the phase transitions and consequently the 
conditions for a self-similar pattern to occur in the magnetization curves. An 
analogous behaviour may also apply for the transport properties of the system 
and this provides a topic for future studies. 
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Figure Captions 

(1) Fig. 1. Magnetization (in units of the saturation magnetization) versus 
applied field H for (a) N = 3, (b) N = 5 and (c) N = 7. The parameters 
are H hl = -1.0 kOe, H ua = 0.5 kOe, H hq = 0.1 kOe and 6 H = 9 ua = 0°. 

(2) Fig. 2. As in Fig. 1 but for H bq = 0.7 kOe. 

(3) Fig. 3. Magnetization (normalized) versus applied field Hq for three dif- 
ferent values of angle $h, taking N = 3, 9 ua = 0° and other parameters 
as in Fie;. 1. 
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